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MAPPING THE CORONARY ARTERIES ON A SPHERE 



The present invention relates to the anatomical imaging arts. It finds 
application in conjunction with the display of the entire coronary tree in the heart context 
and will be described with particular reference thereto. Of course, the technique is also 
applicable to mapping other anatomical regions, such as the cerebral cortex. Although 
5 described with particular reference to CT tomography, it will further be appreciated that the 
invention is equally applicable to other diagnostic imaging techniques which generate two 
and three-dimensional digital diagnostic images including coronary arteries for analysis. 

A variety of imaging modalities can be used to determine structural, 
functional and perfusion characteristics of the heart and coronary arteries. Substantial 
10 research and development have been concentrated on new technologies and methods that 
can aid the diagnosis and analysis of the coronary arteries, and on planning appropriate 
treatment. Much of this effort has been focused on the improvement of the visualization of 
the coronary arteries and the development of computer-assisted analysis of the resulting 
images. 

15 Imaging and analysis of the coronary arteries serve several purposes: 

visualization of the coronary tree, detection and quantification of stenosis, quantification of 
the coronary tree reserve, and analysis of the vessel walls. Traditionally, x-ray angiography 
has been used to detect and quantify stenosis in the coronary arteries. The development of 
multi-slice CT scanner technology having increasingly improved resolution particularly in 

20 the slice direction made the CT imaging of vascular systems attractive for clinical 
applications. Multislice detectors allow several slices to be acquired and processed 
simultaneously enabling more accurate quantification of stenosis. 

Many visualization techniques and quantification strategies are currently in 
use. One approach is the construction of multi-planar reformatted images (MPR) and/or 

25 slab maximum intensity projections (slab-MIP) in optimally chosen planes, so that the 
major segments of the coronary tree are visualized. In addition, the CT's ability of 
acquiring 3D volumetric images can be used to acquire volume rendered images for 
visualizing an isolated coronary arteries tree. However, the current methods for coronary 
arteries inspection, e.g. curved MPR visualization of a single vessel, the planar slab MIP 
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visualization and the volume-rendered visualization of an isolated tree either lack the 
completeness or the context of the coronary tree. 

There is a need for a visualization technique that will allow the volumetric 
display of the entire coronary arteries tree in context, i.e., location, connectivity and 
5 surroundings. The present invention provides a new imaging apparatus and method which 
overcomes the above-referenced problems and others. 

In accordance with one aspect of the present invention a diagnostic imaging 
10 system for displaying a vessel tree is disclosed. A means defines a base surface. A means 
grids the base surface to define pixels. A means projects along a normal of each pixel. A 
means assigns each pixel a grayscale value based on grayscale value of voxels intersected 
by a corresponding normal. 

In accordance with another aspect of the present invention a method of 
15 displaying the coronary arteries tree is disclosed. A base surface is defined. The base 
surface is gridded to define pixels. A normal of each pixel is projected along. Each pixel is 
assigned a grayscale value based on grayscale value of voxels an associated normal 
intersected. A true surface is determined. 

One advantage of the present invention resides in displaying the entire 
20 coronary arteries tree in its context. 

Another advantage resides in using a closed non-planar surface as a base for 
displaying the data of interest. 

Another advantage resides in minimization of distortions in the 
visualization modes. 

25 Another advantage resides in enabling the user to explore the entire 

structure of the cardiovascular system at once. 

Another advantage resides in pennitting the user to adjust the thickness of 
the slab MBP to include more or less data in the image presented. 

Still further advantages and benefits of the present invention will become 
30 apparent to those of ordinary skill in the art upon reading and understanding the following 
detailed description of the preferred embodiments. 
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The invention may take form in various components and arrangements of 
components, and in various steps and arrangements of steps. The drawings are only for 
5 purposes of illustrating the preferred embodiments and are not be construed as limiting the 
invention. 

FIGURE 1 is a diagrammatic illustration of a diagnostic imaging system in 
accordance with the present invention; 

FIGURE 2 is a diagrammatic illustration of a portion of a diagnostic 
10 imaging system in accordance with the present invention; 

FIGURE 3 is an illustration of a base sphere gridded in accordance with the 
present invention and arteries in the vicinity of the sphere; 

FIGURE 4 is an illustration of a coronary arteries tree in a true amorphous 

surface form; 

15 FIGURE 5 is an illustration of a coronary arteries tree in a globe surface 

form; and 

FIGURE 6 is a 2D map of a coronary arteries tree. 



20 With reference to FIGURE 1, a diagnostic imaging apparatus 10 generates 

electronic diagnostic image representations. In the preferred embodiment, the diagnostic 
imaging apparatus 10 includes a CT scanner which generates a three-dimensional 
volumetric image representation which is made up of voxels whose position is defined 
along orthogonal axes. Other diagnostic imaging systems, such as magnetic resonance 

25 imaging, PET imaging, SPECT imaging, and other diagnostic techniques which generate 
analogous three-dimensional diagnostic images are also contemplated. 

In the preferred embodiment, the CT scanner includes a stationary gantry 
12, in which a rotating gantry 14 is mounted. The rotating gantry carries an x-ray tube 16 
and a two-dimensional array of detectors 18, which are diametrically opposed to each other 

30 across a scan circle or an examination region 20. A subject support 22 supports a region of 
interest of the subject in the examination region 20. Prior to imaging, the subject is 
preferably injected with a known contrast agent that produces the voxels of the known 
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intensity in the vessels, preferably, at the intensity higher than that of tissue. A 
longitudinal drive motor 24 moves the subject support 22 longitudinally through the 
examination region 20. Preferably, the subject support 22 longitudinally reciprocates the 
subject as the rotating gantry 14 rotates continuously for spiral scanning. The imaging 
5 apparatus generates a volumetric image representation of a transverse volume of the 
subject, which includes the subject's heart, or other region of interest. Alternately, the 
subject support 22 can be stepped and data can be collected along a series of parallel, 
transverse slices. Although the detector array 18 is illustrated as rotating with the rotating 
gantry 14, the detector array 18 may optionally be mounted as a continuous ring on the 

1 0 stationary gantry 1 2 . 

A three-dimensional reconstruction processor 26 reconstructs the output 
signals from the detector array 18 in accordance with the angular position of each detector 
element, the angular position of the x-ray tube 16, and the longitudinal position of the 
subject support 22 at the time of sampling into one or more three-dimensional image 

15 representations. The reconstructed image representation is stored in a volumetric image 
memory 28. 

With continuing reference to FIGURE 1 and further reference to FIGURES 
2 and 3, a base surface means or processor 30 approximates a base surface 32, which is 
preferably a sphere. Of course, other surfaces such as ellipsoids, anatomical heart 

20 approximations, and the like, are also contemplated. The approximation is based on the 
assumption that the coronary arteries lie on a substantially smooth closed surface. More 
specifically, a volume selecting means 34 selects a volume of data from the volume 
memory 28 which corresponds to the region of interest, e.g. heart, brain, or other organ. A 
centerlines determining means 36 finds the centerlines 38 of the vessels in the selected 

25 volume by one of known techniques. A best fitting surface means or process 40 draws a 
best fitted surface, preferably the sphere or ellipsoid, to the determined centerlines. Some 
of vessels will lie above the surface of the base surface 32, some of the vessels will he 
underneath it, and some will have a portion above and a portion below the surface. The 
sphere or ellipsoid is rotated such that the axis of rotation is substantially parallel to a long 

30 axis of the left ventricle. In performing the best fitting process 40, a centerlines coordinates 
converting means 42 converts the centerlines coordinates to spherical coordinates 
according to equations 1. 
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<p = a tan [Z / V (X 2 + Y 2 )] Equations 1 

A = atan[Y/X] 

A = [V(X 2 + Y 2 )/cosp]-R 

5 

where <p is the latitude; X is the longitude; h is the distance from the sphere; X, Y, Z are the 
Cartesian coordinates of a centerline point; and R is the radius of the sphere. 

With continuing reference to FIGURES 1, 2 and 3, a gridding means 50 
spreads a grid over the base surface 32 to grid the sphere, or other determined surface, into 

10 pixels 52. The base surface 32 is stretched or shrunk along the sphere normals to fit the 
true form of the vessels wherever a centerline is located such that the stretched surface lies 
on the centerlines. With further reference to FIGURE 4, the stretching and shrinking of the 
sphere to fit the centerlines is performed by a true surface detenriioing means 54 which 
builds a true surface 56 representing a 'True Form" undistorted mode visualization, in 

1 5 which the coronary arteries tree in its context is revealed over a true form surface running 
through its vessels. More specifically, a surface adjusting means, 58 receives the 
centerlines information from the centerline means 36 and moves the location of the 
spherical base surface 32 along the normals, perpendicularly to the surface, to the point 
where the normals intersected associated centerlines of the vessels. The remainder of the 

20 spherical base surface is mapped in and out to define the smoothly varying true form 
surface. The resulting true surface is the arbitrary, amorphous surface, as seen in FIGURE 
4. 

A thickness determining means 60 determines a wall thickness of the true 
surface 56. Preferably, a user supplies predetermined thickness value to the thickness 

25 determining means 60, e.g. the thickness of the hollow true form volume. Alternatively, 
the thickness is determined by the thickness determining means 60 based on specified 
criteria, such as vessel diameter, or other characteristics of the subject. A normal 
projecting or tracking means 62 projects a normal from each pixel 52 orthogonally with 
respect to the true surface 56 in both directions, e.g. outward and inward from the true 

30 surface 56, to the limits of the wall thickness determined by the thickness determining 
means 60. A grayscale assigning means 70 examines the intensity of each voxel along the 
normal and assigns a grayscale or intensity value to the corresponding pixel of the true 
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surface 56 to be displayed in the resultant image. A screening means 72 screens the 
grayscale values of voxels, intersected by each normal, based on predetermined criteria to 
select a grayscale intensity value which will be displayed for the corresponding pixel. In 
the preferred embodiment, the screening means 72 selects a maximum intensity value 
5 along each normal and stores it in a maximum intensity image pixel memory 74. The 
cycle repeats until all the normals are screened and a grayscale value of each pixel is 
determined. 

With continuing reference to FIGURE 1 and further reference to FIGURE 
5, a globe means 80 presents a "Globe" mode of visualization, in which the coronary 

10 arteries tree is shown on a sphere which is best fitted to the determined amorphous true 
surface 56. More specifically, a draping means 82 drapes or projects the grayscale values, 
which have been assigned by the grayscale assigning means 70 to the true form surface 
onto the spherical base surface by projecting the voxels along the normal to the sphere 
corresponding to each spherical or globe surface pixel. The resulting surface is a globe like 

15 structure 84 and may be rotatably examined as a globe. 

To display the "True Form" or "Globe" on a conventional monitor 90 
having a two-dimensional display, a video processor 92 projects a normal from each pixel 
of the display. Each display pixel is assigned the grayscale value of the true surface or the 
globe surface that it intersects. It will be noted that surface pixels closest to the viewing 

20 screen are displayed substantially undistorted, while surface pixels toward the poles and 
the periphery become progressively more compressed to give the viewer the impression of 
viewing a 3D object on the 2D viewing screen. Optionally, the 3D effect can be enhanced 
by assigning an illumination direction and adding shading. 

With continuing reference to FIGURE 1 and further reference to FIGURE 

25 6, a 2D means 100 maps the globe 84 or a selected portion of it on a flat surface analogous 
to the way that the round surface of the earth is projected onto flat maps. More specifically, 
a 2D gridding means 102 spreads the original grid of the base surface 32 on a 2D x-y 
surface. Typically, equatorial surface pixels are substantially undistorted, while more polar 
surface pixels are stretched to maintain spatial relationships. A matching means 104 

30 matches each coordinate of the base surface 32 to corresponding coordinates on the x-y 
viewing plane. A 2D grayscale processor or means 106 assigns each pixel on the x-y 
viewing surface a grayscale value of the corresponding pixel from the globe surface 84. 
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A mode selection means 108 receives an input from a user and selects a 
visualization mode for processing. The selected visualized data is processed by the video 
processor 92 and displayed then on a monitor 90, A user input device 110, such as a 
keyboard, enables the user to control the video processor 92 to display selected projections, 
5 rotate the globe or true surface, and input necessary information as was discussed above. . 

In one alternate embodiment, the true form surface and its wall thickness 
are defined, but then data is projected directly onto the spherical (or elliptical) base surface. 
In another alternate embodiment, the two dimensional image is mapped directly from the 
true form image rather than the globe image, which was derived from the true form image. 
10 The invention has been described with reference to the preferred 

embodiments. Modifications and alterations will occur to others upon a reading and 
understanding of the preceding detailed description. It is intended that the invention be 
construed as including all such modifications and alterations insofar as they come within 
the scope of the appended claims or the equivalents thereof. 
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